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The optical conductivity ol LaFeAsO, BaFe2As2, SrFe2As2, and EuFe2As2 in the spin-density 
wave (SDW) state is investigated within density functional theory (DFT) in the framework of spin- 
polarized generalized gradient approximation (GGA) and GGA-I-U. We find a strong dependence 
of the optical features on the Fe magnetic moments. In order to recover the small Fe magnetic 
moments observed experimentally, GGA-t-C/cff with a suitable choice of negative on-site interaction 
Ueff = U — J was considered. Such an approach may be justified in terms of an overscreening 
which induces a relatively small U compared to the Hund's rule coupling J, as well as a strong 
Holstein-like electron-phonon interaction. Moreover, reminiscent of the fact that GGA-t-f/off with a 
positive f/eff is a simple approximation for reproducing a gap with correct amplitude in correlated 
insulators, a negative UeS can also be understood as a way to suppress magnetism and mimic 
the effects of quantum fluctuations ignored in DFT calculations. With these considerations, the 
resulting optical spectra reproduce the SDW gap and a number of experimentally observed features 
related to the antiferromagnetic order. We find electronic contributions to excitations that so far 
have been attributed to purely phononic modes. Also, an orbital resolved analysis of the optical 
conductivity reveals significant contributions from all Fe 3d orbitals. Finally, we observe that there 
is an important renormalization of kinetic energy in these SDW metals, implying that the effects of 
correlations cannot be neglected. 

PACS numbers: 74.70.Xa. 71. 15. Mb. 74.25. Gz, 74.25.Jb 



I. INTRODUCTION 



In several families of the iron pnictides, high temper- 
ature superconductivity (SC) emerges in close proximity 
of an antiferromagnetic (AF) ground state with stripe- 
type order. ^ The AF transition is induced by a spin- 
density wave (SDW) instability below a critical tempera- 
ture TsdWj which is either preceded or coincidental with 
a structural transition from a tetragonal phase to an or- 
thorhombic phasci^ Upon doping^ or application of pres- 
sure^, the AF order is suppressed and superconductivity 
emerges. Like lattice vibrations, AF fluctuations are ca- 
pable of producing an attractive interaction necessary for 
the creation of Cooper pairs. It is therefore widely be- 
lieved that the AF fluctuations drive the superconducting 
instability in these systems. As for the nature of the AF 
state, the iron pnictides show signatures of both electron 
itinerancy and local magnetism. Some experimental^^— 
and theoretical works based on density functional theory 
(DFT)i^-i^ and dynamical mean-field theory (DMFT)!^ 
as well as a combination of these two methods^i^i favor 
an itinerant scenario with rather moderate correlation 
strength, as reflected by the metallic nature of the com- 
pounds. On the other hand, other authors point out the 
effects of strong correlations^"— like the renormaliza- 
tion of the kinetic energy of the electrons j^^S 

Many features of the electronic structure of the iron 
pnictides are directly reflected in their optical proper- 
ties: the low-frequency region of the conductivity spec- 
trum is governed by the itinerant carrier contribution and 
directly shows the effect of correlations in the area under 



the Drude region which is proportional to the electron's 
kinetic energy; the infrared regime above the Drude peak 
is dominated by gap features induced by either the SDW 
gap or the superconducting gap; finally, the visible part 
of the spectrum reflects the band structure in the normal 
state. Consequently, a number of experimental studies 
have been performed on the optical properties of the iron 
pnictides, in the normal state, the SDW state, as well as 
in the SC state,^iii'22,-^i^ The SDW state is character- 
ized by (i) the appearance of a peak in the optical con- 
ductivity at the SDW gap frequency, (ii) an anisotropic 
dc response and (iii) an almost isotropic response in the 
infrared and optical region of the spectrum. On the other 
hand, while several theoretical works based on local den- 
sity approximation (LDA)?^? and LDA-hDMFT^^'^" have 
been done on the paramagnetic phase, there is still a 
lack of DFT calculations for the optical conductivity in 
the SDW state. 

In this work we report optical studies in the framework 
of density functional theory on four iron pnictides in the 
SDW state, namely the 1111 compound LaFeAsO, and 
the 122 compounds ^£;Fe2As2 (^£;=Ba,Sr,Eu). This 
provides an insight into the microscopic origin of the op- 
tical features in the SDW state and allows for an as- 
sessment of DFT regarding its applicability to the iron 
pnictides. 



II. COMPUTATIONAL DETAILS 

We performed electronic structure calculations with 
the full potential linearized augmented plane wave 



(FLAPW) method as implemented in WIEN2K.35 The self- 
consistency cycle employed 2048 k points in the full Bril- 
louin zone (FBZ) using the generalized gradient approx- 
imation (GGA) in the Perdew-Burke-Ernzerhof variant 
for the exchange correlation potential;"^^ the optical prop- 
erties were evaluated with 16384 k points in the FBZ, 
the number of k points in the irreducible Brillouin zone 
depends on the symmetries of the respective material. 
Experimental lattice parameters and atomic positions 
were used, from Refs. [STl - liol for LaFeAsO, BaFe2As2, 
SrFe2As2, and EuFe2As2, respectively. All calculations 
were performed in the scalar relativistic approximation. 
For the optical properties, the optics'*^ code package in 
WIEN2K was modified to allow for an orbital character 
resolved analysis. 

We are working in the framework of 'GGA-I-U' where 
'[/'(= Ucff = U — J) describes the competition between 
the (spherically averaged) on-site Coulomb interaction U 
and the (spherically averaged) on-site exchange coupling 
J (within the Fe 3d subshell for the iron pnictides). In 
this context, the atomic limit double counting correction 
according to Refs. |43 andlsil was applied. With this dou- 
ble counting correction, the expression for the correction 
to the GGA functional for the orbital m reads 



U ~J 
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'^m<7 \^ 
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where rima is the spin-projected occupation in orbital m. 

In order to allow for the stripe-type AF order, we con- 
sider a doubled (-^72 x y/2 x 1) unit cell with AF order 
along the a axis of the supercell and ferromagnetic (FM) 
arrangement along the b axis {i.e. the supercell is rotated 
45° with respect to the original unit cell) , as observed ex- 
perimentally. In the following, the orbital characters and 
dielectric tensor components are labelled with respect to 
the coordinate system of this supercell. Spin-polarized 
calculations with AF order are labelled with 'GGA(AF)' 
('GGA-I-U(AF)', respectively). 

The linear optical response of the electronic systeroSE 
is calculated in the random phase approximation (RPA) , 
with the Kohn-Sham orbitals mimicking the bare elec- 
tronic states of the electron-hole excitations in the RPA 
formalism. In this framework, the complex- valued dielec- 
tric function e is calculated as 
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where Q.c is the volume of the unit cell, q the wavevec- 
tor of the incoming light, /o the Fermi distribution, and 
{n, fc} labels the Kohn-Sham (KS) orbital in band n with 
crystal momentum k. The transition matrix element is 
given by Mn^n'{k,q) = {n,k\e^'^'^'^\n' ,k + q). As ex- 
pressed by limg_^o, only direct transitions fc/ — ki are 
considered here since even for high energies, the wave- 
length of the light is large compared to the lattice di- 
mensions. In this limit, perturbation theory for k p can 



be applied (p being the momentum operator) in order to 
obtain expressions for \n,k + q) and En^k+q in terms of 
|n, fc) and En^k within the dipole approximation. This 
yields for the dielectric tensor— 
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where the second and third term describe the intraband 
and interband contributions, respectively. Since T = 
is considered, the intraband contribution is restricted to 
states at the Fermi energy Ep, whereas the band indices 
v{c) run over all occupied (empty) states, pi is the ma- 
trix element of the momentum operator along the elec- 
tric field polarization of the incoming light, Pi-^n,n',k = 
{n,k\pi\n' ,k). Spin-orbit coupling is not taken into ac- 
count in our calculations, and therefore all off-diagonal 
elements i ^ j vanish due to the orthorhombic symmetry 
of the crystals. 

Turning the sum in Eq. ^ into an integral over the 
first Brillouin zone and using the Dirac representation 
■^TTj- — V- ^i7TS{x), the imaginary part of the interband 
contributions to the dielectric tensor components can be 
expressed as 



y^ / \Pi;c,v,k\'^S{Ec,k-Ey^k-fl^). 
„ ^ Jk 
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(4) 
In this work, we focus on the analysis of the real part 
of the optical conductivity. 



Rea.ij{uj) = -—Imeij{uj). 
47r 



(5) 



For metals, the intraband contribution has the well- 
known Drude shape 
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r^2 
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where the plasma frequency ujp determines the spec- 
tral weight of the Drude peak via Jp dui au^uj) — ^, 
whereas its width is given by the carrier scattering rate 
r — l/r with the lifetime r. 

The presence of the Drude peak depends on the metal- 
licity of the material under consideration and in general 
needs to be confirmed by an explicit calculation of the in- 
traband transitions. However, as only direct transitions 
are considered here, the intraband part of Eq. ([3]) is sin- 
gular at w = and vanishes at finite frequencies, both 
for the dielectric function as well as for the conductivity. 
Instead, the squared plasma frequency for band n and 
spin a is defined as 

ujl{n,<j)^ / \{n,k,(j\p\n,k,a)\^5{En,k.a - Ey) 

(7) 



with the total squared plasma frequency uj^ = 
X]„ o.'^p('^,cr), and the Drude peak is considered as an 
approximation to the intraband conductivity, i. e. we con- 
sider only one overall Drude peak, despite the possibil- 
ity of two or more Drude peaks due to multiple bands 
crossing the Fermi surface. For the iron pnictides, this is 
justified as they show metallic behaviour also below the 
SDW transition. Consequently, only the carrier scatter- 
ing rate is left as free parameter which cannot be de- 
termined within a DFT approach. For our purposes, 
r is chosen close to experimental values^ in the SDW 
state. Given these considerations, all relevant informa- 
tion about the Drude peaks is contained in the plasma 
frequency uip. 



III. RESULTS 
A. Density of States 

Preliminary to the analysis of the optical properties of 
the iron pnictides, we first start with the discussion of the 
density of states (DOS) since the dielectric function is es- 
sentially the joint density of states (jDOS) weighted with 
the transition matrix element (see Eq. (|4])). In particular, 
we find that the antiferromagnetic order and the magni- 
tude of the magnetic moments crucially affect the optical 
properties which is reflected by the (partial) opening of 
a gap in the DOS. 

For the iron pnictides, DFT calculations are known to 
strongly overestimate the magnetic moments on the Fe 
atoms in the SDW state (« 2.0/^b with GGA(AF) com- 
pared to 0.4/XB^-0.6/iB^ for LaFeAsO). As wiU be dis- 
cussed in the following, this results in a too large SDW 
gap which in turn shifts the SDW features in the optical 
properties to erroneously high energies. In order to repro- 
duce correctly the SDW gap, we employ GGA+U with 
Ucs negativ o^'^^ to suppress the overestimated magnetic 
moment from GGA(AF). On the one hand, the nega- 
tive C/off can be viewed as a result of an overscreening of 
the on-site Coulomb interaction which leads to a smaller 
value of U compared to the on-site exchange interaction J 
or a strong Holstein-like electron-phonon coupling which 
gives a frequency-dependent negative contribution to the 
on-site Coulomb interaction. On the other hand, similar 
to the application of positive U at the mean-field level 
in GGA-I-U calculations to reproduce the correlated gap 
in Mott insulators, the negative C/cff can be understood 
as a simple way to suppress the tendency to magnetism 
within the mean-field approximation and mimic the ef- 
fect of quantum fluctuations ignored in DFT. This can 
be read off from the correction to the GGA functional in 
Eq. ^ which energetically favors Uma — 1/2 at negative 
UcS, i.e. the paramagnetic case, and conversely penalizes 
(spin-projected) integer occupations, i.e. the fully polar- 
ized states. Note that the exact form of the GGA-I-U 
functional and thus the preference of low-spin or high- 
spin states depends on the double counting correction. ^^^ 



20- SrFe/s^ 




FIG. 1. (Color online) DOS of SrFe2As2 for different val- 
ues of C/cff. oj = corresponds to the Fermi energy. As is 
well known, the Fe 3d subshell dominates the DOS around 
the Fermi energy. The contributions from Sr and As, and 
the other Fe orbitals are small. The magnetic moments from 
top to bottom are m = 0, m — 1.98 fiB, rn = 1.30 /iB, and 
m — 0.43 ^B, respectively. See text for a discussion of the 
transitions labelled by ti and t2. 



A comparison of the role of U and J with their model 
counterparts is therefore not generally valid but needs to 
take into account the double counting correction under 
consideration. Also, the magnetic moment needs to be 
further reduced below its experimental value due to the 
RPA employed in our calculations of optical conductivi- 
ties, as discussed later. Therefore, the negative C/off used 
in this paper is not comparable to realistic values of U 
and J for these systems due to the approximations we 
employed here. 

Fig. [T] shows the evolution of the density of states for 
different values of Ucs in SrFe2As2. Compared to the 
non-spin-polarized GGA DOS, a suppression of the DOS 
around Ep is clearly visible in the antiferromagnetic cal- 
culations. For the parameter regimes shown in Fig. [TJ 
only partial opening of the gap is observed as some DOS 
persists around E-p; full opening of the gap occurs at posi- 
tive values of Ues, i.e. for sufficiently strong Hubbard U. 
In contrast, for pronounced negative C/off the magnetic 
moment is reduced and the SDW gap narrows accord- 
ingly. For C/off = — 1.9 eV and m = 0.43 /iB, the energy 




:(b) 


BaFe^As^ 






- 


A- 


! 


J\ 


■ 


/K\ 


ir 


/ v — 


T 


/ V 




X > 



0.1 0.2 

to(eV) 



0.3 



(f) 


^^xz/xz 


1 l/Z/X|/ 


t^ 


\ ^Mxy/xz 


/ 


^ 1 l/Z/KZ 


1 1^ 


\ \ ^xy/yz 


/// 


V^ \ ^M^^/^y 


§ 





0.2 0.25 

(0(eV) 



0.1 



0.15 
m(eV) 



0.2 




5 


[ 


(d) LaFeAsO 


4 




A 


"oT 


■l 


aJ><)l. 


=o3 




/Tv — ^ 


e 2 


IP 


K/ V 


1 

n 





0.1 0.2 

to(eV) 



0.3 



I \vz/xv 
I \xz/vz 



0.15 0.2 

(0(eV) 



0.15 0.2 

co(eV) 



FIG. 2. (Color online) Real part of the optical conductivity of different iron pnictides in the SDW state, for UeS = —1.9 eV: 
(a),(e) SrFe2As2 (experimental data taken from Ref. ,5), (b),(f) BaFe2As2 (Ref. |^), (c),(g) EuFe2As2 (Ref. [sj), and (d),(h) 
LaFeAsO (Ref. |11| ). No (empirical) Lorentz broadening has been applied; in order to facilitate the comparison with the 
damped/broadened experimental data, the experimental results have been rescaled by a factor of two. The orbital character 
contributions are sorted by their contribution to the total conductivity, integrated over the energy range shown in the re- 
spective plot, with the largest contribution at the bottom (for clarity, only the first seven contributions are displayed; smaller 
contributions have been omitted). 



gap is A w 0.17 eV. Almost over the whole energy range 
shown in the figure, the DOS is strongly distorted and 
shifted by the inclusion of the antiferroniagnetic order 
compared to the non-spin-polarized GGA DOS. How- 
ever, the strong suppression of the magnetic moments 
for [/off = —1.9 eV renders the DOS almost equal to the 
GGA DOS again, except in the region around the Fermi 
energy. Thus, the optical properties can be expected to 
be close to the GGA case in the higher energy regions of 
the spectrum. 



B. Optical conductivity 

As for the in-plane optical conductivity of the iron 
pnictides in the SDW state, experimental investigations 
on single crystals^^i^i^ consistently report a number 
of common features, although the exact location of the 
features depends on the material under investigation. 
These features are: metallic behaviour, i.e. the presence 
of a Drude-like conductivity at low frequencies (< 100 
cm~^), a sharp peak at the SDW gap frequency (« 1000 
cm~^-1500 cm~^), and a broad, less pronounced peak in 
the midinfrared region (w 5000 cm~^~6000 cm~^) which 



almost doesn't depend on temperature and is present 
also above Tsdw- The peaks are associated with a sup- 
pression of the spectral weight at lower energies (below 
« 1000 cm-i for the SDW peak, below « 5000 cm^^ for 
the high-energy peak) which leads to a spectral weight 
transfer from lower to higher energies. Note that in the 
normal state, the tetragonal symmetry leaves only two in- 
dependent components in the dielectric tensor and thus 
the conductivity is defined as aaa = o'bb for the Fe in- 
plane directions a and b, and ace for the out-of-plane di- 
rection c perpendicular to a and b. In the SDW state, the 
stripe-like AF order introduces an anisotropy which lifts 
the degeneracy of the two in-plane components. How- 
ever, no substantial anisotropy is found experimentally. 

In Fig. [2] we present an overview of the in-plane opti- 
cal conductivity in the low-frequency region as obtained 
by our calculations. Figs. [2] (a)-(d) show a comparison 
of the experimental results with the calculation methods 
introduced in the previous section; the GGA-t-U(AF) re- 
sults exhibit the SDW peak as the most prominent fea- 
ture, located at the SDW gap frequency. For all com- 
pounds, the SDW peak emerges at the experimentally 
determined frequency only for a suitable negative value 
of Ucff, and otherwise moves to higher frequencies (see 




FIG. 3. (Color online) Optical conductivity in the higher en- 
ergy region for (a) LaFeAsO, (b) BaFe2As2, and (c) SrFe2As2. 
Same experimental references and values for UcS as in Fig. (2] 
The arrows indicate the position of the respective SDW peak. 



Fig. Hfor SrFeaAsa). The value U^s = -1-9 eV was 
approximately determined by denianding a correct SDW 
peak position. Fig. [3] displays the optical conductivity 
for LaFeAsO, BaFe2As2, and SrFe2As2 over a larger en- 
ergy range and includes some damping. As can be seen 
there, the GGA results vifithout SDW order don't show 
any significant peak in the optical conductivity up to 
« 1 eV. This can be readily read off from the density of 
states, where the GGA DOS is essentially depleted up to 
1 eV above Ep. At higher energies, there is an obvious 
disagreement between the DFT results and the experi- 
mental data, for both GGA-l-U(AF) and GGA. Whereas 
we already mentioned the midinfrared peak around 0.6 
eV, the experimental data basically doesn't show any 
structure above 1 eV. This could be described if quan- 
tum fluctuations are treated properly and consequently 
a frequency-dependent self-energy is involved. However, 
as the focus of this work is the analysis of the features in- 
duced by the antiferromagnetic ordering, we concentrate 
on the low energy features as indicated by the arrows. 

As for the anisotropy due to the stripe-type SDW 
symmetry-breaking, we find that it is strongly present 
in the GGA(AF) calculation {i.e. for UcS — 0) and for 
moderately negative values UoB ^ — 1-5 eV, where atb - 
the conductivity along the FM axis - shows a peak at 
significantly higher frequencies than aaa- In fact, this 
anisotropy can be mapped to the DOS in Fig. [U where 
ti denotes the transition(s) which dominate the peak in 
abb and t2 the ones dominating aaa- As Ucs is decreased, 
the anisotropy is gradually suppressed and for the regime 
shown in Fig. [21 almost no anisotropy is present any- 
more in SrFe2As2, BaFe2As2 and EuFe2As2, in accor- 
dance with the experimental observations. This can be 
understood by the fact that while dxz and dyz are degen- 
erate in the high temperature tetragonal phase, below 



the SDW transition temperature, different band split- 
tings and band shiftings are produced by different mag- 
netization and occupation number on these two orbitals. 
Therefore, the excitations along x and y direction be- 
come inequivalent. As UcS is decreased, the magnetiza- 
tion and consequently the symmetry-breaking between 
dxz and dyz orbitals are suppressed, leading to smaller 
differences between dxz and dyz orbitals and accordingly 
in the excitations along x and y direction. An anisotropy 
can still be seen in the weak double peak structure of 
LaFeAsO and SrFe2As2 in Figs. [11(d) and (h), respec- 
tively (a) and (e) (the conductivity shown in the plots 
is given by cr = {aaa + crbb)/"^)'- e.g. in LaFeAsO, the 
dxy contribution is of roughly equal size in both sub- 
peaks, whereas the dxz contribution is mainly present in 
the first sub-peak and the dyz contribution in the second 
one. Consequently, the second sub-peak is dominated by 
abb, J-e. by the conductivity along the FM direction. We 
want to point out that also this sub-peak is directly re- 
lated to the AF order as no peak is present in a purely 
ferromagnetic calculation (with the same absolute value 
of the magnetic moments). 

The observed midinfrared peak at saO.6 eV, which is 
visible in Fig. [3l for LaFeAsO and BaFe2As2, respec- 
tively, is absent in all calculations, including the non- 
spin-polarized GGA calculations. This is in contrast to 
the experimental data and needs to be investigated fur- 
ther. 

Surprisingly, the calculated optical spectrum of 
SrFe2As2 shows an experimentally observed low- 
frequency feature at 31 meV (see Fig. HJa)) that has 
so far been attributed to phononic excitations: while 
there is no discussion in the respective experimental work 
Ref. [^ for SrFe2As2, for EuFe2As2 it is argued in Ref. [S^ 
that the phononic mode - there at 32 meV - exhibits 
a coupling to the electronic background only for high 
temperatures while for low temperatures the Lorentz line 
shapes indicates a purely phononic, non-interacting exci- 
tation. However, our results show that this peak can also 
be obtained from the electronic band structure alone in 
SrFesAss. 

In Figs. [51(e)-(h), the decomposition of the SDW peak 
into Fe 'id orbital characters of the initial and final states 
is shown. The orbital resolved optical conductivity reads 
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where A™^. formally denotes the relative weight of the 
Fe orbitals with magnetic quantum number m in the KS 
orbital |n, fc). Note that the momentum operators only 
couple states with different parity, reflected by the well 
known selection rule ^l = lf — U ~ ±.1 for dipolar transi- 
tions. Dipolar transitions among Fe Sd states are there- 
fore forbidden in the atomic limit. For the iron pnictides, 
the states around the Fermi energy are Fe id dominated 
with some As Ap contribution. This hybridisation al- 
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FIG. 4. Dependence of the SDW peak position in the optical 
conductivity and the magnetic moment on the value of UcS 
for SrFe2As2. The thin lines indicate experimental values: 



Mcxp = 1.01/iE 



SDWc 



0.18 eV-^ 



lows for transitions between initial and final states which 
are both Fe 3d dominated but owe their finite transition 
strength to Fe 3d O As ip transitions. This is taken into 
account in Eq. ([5]) where the full wavefunction |n, fe) is 
used for the calculation of the matrix element but the 
resulting dielectric function is projected on the Fe 3d 
subspace. This projection in particular omits the con- 
tribution from the interstitial region which has a consid- 
erable DOS in the considered energy range but cannot be 
assigned an orbital character. However, as can be seen 
from the comparison of the SDW peaks in Figs. [2] (a)- 
(d) - which show the total optical conductivity - to the 
close-up in (e)-(h) - which only show the Fe 3d contribu- 
tion -, the Fe 3d part resembles very well the structure of 
the total conductivity. Therefore the contributions from 
the interstitial region as well as from the other atoms are 
neglected for the orbital character analysis. 

In the iron pnictides, the Fermi surface is crossed by 
multiple bands. Consequently, also the optical properties 
have multiband nature and one can expect contributions 
from several orbital characters. This is confirmed by our 
calculations where we observe no overly dominating char- 
acter component in any part of the spectrum. The SDW 
peak structure differs from the rest of the spectrum in 
that it mainly contains t2g character components; all t2g 
components are larger than any Cg component. 

Finally, we analyze the low-energy Drude region of 
the spectrum as characterized by the plasma frequency. 
Since the SDW gap opening is only partial, the Drude 
peak is still present - expressed by a finite plasma fre- 
quency ~ both experimentally and in our calculations. 
The ratio of the kinetic energies which equals the ratio 
of the squares of the plasma frequencies, iiTexp/i^band — 
(a;p''P)^/(a;p™'^)^, is commonly taken as a measure for 
the renormalization effect from the electronic correla- 
tions compared to band structure calculations. As given 
in Table m the non-spin-polarized GGA value for, e.g., 



TABLE I. Optical parameters of the investigated compounds 
as obtained from the different calculation methods, cjp is 
the plasma frequency in a(6)-direction. 



-~^ Compound Calc. 



[eV] io^ [eV] # bands m [^b] 
at Ef 



LaFeAsO 



BaFe2As2 



SrFe2As2 



EuFe2As2 



GGA 

GGA(AF) 

GGA+U(AF), 
(7eff = -1.9eV 

GGA 

GGA(AF) 

GGA-hU(AF), 
[/off = -1.9eV 

GGA 

GGA(AF) 

GGA-HU(AF) 
[/eff = -1.9eV 

GGA 

GGA(AF) 

GGA-hU(AF) 
f/eff = -1.9eV 



2.25 
0.87 
1.74 

2.62 
0.62 
1.66 

2.79 
0.64 
1.50 

2.96 
0.84 
1.95 



2.25 5 
0.82 2 
1.17 2 



2.62 
0.65 
1.53 

2.79 
0.64 
1.27 



2.96 5 
0.99 2 
2.04 3 





1.98 

0.29 



1.98 

0.30 



1.98 

0.43 



1.76 

0.34 



SrFe2As2 is u!p « 2.79 eV, whereas the experimental 
value in the normal state is uip sa 1.7 eV at 300 K^ 
yielding K^xp/Kband ~ 0.37. In the SDW state, the 
experimental plasma frequency is strongly reduced due 
to the removal of itinerant carriers from the Fermi sur- 
face by the opening of the SDW gap, to ujp ~ 0.59 eV 
at 10 K.^ Likewise, in our calculations, the inclusion of 
the SDW order reduces the number of bands crossing 
the Fermi surface, which also significantly reduces the 
plasma frequencies. The renormalization due to the cor- 
relations still persists, though: with A'cxp/^band ~ 0.15 
at Ucs = — 1.9 eV it is even more pronounced in the SDW 
state. 

In Table U as well as Fig. |31 one also notices that 
with our choices of C/cff in the GGA+U(AF) calcula- 
tions, experimental SDW peak positions are well repro- 
duced, while the calculated magnetic moments are all 
smaller than the corresponding experimental values, for 
example, the calculated magnetic moment for SrFe2As2 is 
0.43 /is while it is l.Ol/is from neutron diffraction exper- 
iments'*^''*^. The discrepancy is due to the fact that the 
employed RPA scheme neglects correlation effects, like in- 
teractions between electrons and holes which may reduce 
the excitation energy and would therefore shift the SDW 
peak towards lower frequencies. In order to account for 
such a shift induced by correlation effects, a smaller value 
of the magnetic moment compared to the experimental 
one has to be used, which similarly shifts the SDW peak 
towards a lower position. For LaFeAsO, the discrepancy 
between the calculated magnetic moment, with which 
the experimentally observed SDW peak position is re- 
produced, and the experimental one is less pronounced, 
indicating less correlation strength in LaFeAsO than in 



the other three compounds which is consistent with the 
DFT downfolding results of Ref. l48|- 

For all compounds except EuFe2As2, the plasma fre- 
quencies and thus the dc conductivities exhibit a notable 
anisotropy between the a and b axis. Interestingly, the 
anisotropy develops differently in the 1111 compound and 
the 122 compounds where the 1111 compound always fea- 
tures a higher conductivity along the a direction. In con- 
trast, (t'' > (j"^ in the 122 compounds for Ucff = which 
is expected because of the larger lattice constant along 
the a axis and the orientation of the SDW vector along 
a. Unexpectedly, this no longer holds for negative Ucff 
where we find cr° > a'' also for the 122 compounds. This 
is in agreement with recent experiments on underdoped 
BaFe2As2 where the magnetic domains were (partially) 
detwinned in a magnetic field in order to reveal the in- 
plane anisotropy in the resistivity.— 



IV. CONCLUSIONS 

In summary we have demonstrated that DFT is capa- 
ble to reproduce a number of features associated with the 
SDW state in the iron pnictides. However, this comes at 
the cost of a negative Ucs in the context of GGA-I-U cal- 
culations. The negative U may be understood as a way of 



simulating a strong screening of the Coulomb interaction 
and a large electron-phonon coupling^"— i^"—. Also it 
can be viewed as a route to mimic the effects of quan- 
tum fluctuations at the mean-field level and therefore 
as a driving force for suppressing the overestimated Fe 
magnetic moments obtained from GGA(AF). Concern- 
ing the size of the magnetic moments and thus the posi- 
tion of the SDW peak, the negative Ucff therefore needs 
to be interpreted as a fitting parameter. However, the 
agreement with the experimentally observed isotropy of 
the in-plane conductivity and the appearance of the low- 
frequency phononic-electronic excitation at exactly the 
right frequency indicate the reliability of the approach. 
Therefore, we conclude that the used LDA-I-U framework 
with negative Ucs mainly reduces the magnetic moment 
but doesn't distort the overall band structure too seri- 
ously. This suggests that this method is surprisingly well 
suited for the description of the SDW state in the iron 
pnictides. 
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